Critical properties of the half-filled Hubbard model in three dimensions 
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By means of the dynamical vertex approximation (DFA) we include spatial correlations on all 
length scales beyond the dynamical mean field theory (DMFT) for the half-filled Hubbard model 
in three dimensions. The most relevant changes due to non-local fluctuations are: (i) a deviation 
from the mean-field critical behavior with the same critical exponents as for the three dimensional 
Heisenberg (anti)-ferromagnet and (ii) a sizable reduction of the Neel temperature (Tjv) by ~ 30% 
for the onset of antiferromagnetic order. Finally, we give a quantitative estimate of the deviation of 
the spectra between DFA and DMFT in different regions of the phase-diagram. 
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Almost 50 years after the invention of the Hubbard 
model [l| and despite modern petaflop supercomputers, a 
precise analysis of the criticality of this most basic model 
for electronic correlations has not been achieved so far, at 
least not in three dimensions. Dynamical mean field the- 
ory (DMFT) [H-Q was a big step forward to calculate the 
three dimensional Hubbard model since the major con- 
tribution of electronic correlations, i.e, the local one, is 
well captured within this theory. Local correlations give 
rise to quasiparticle rcnormalization, the Mott-Hubbard 
transition, magnetism, and even more subtle issues such 
as kinks in purely electronic models However, non- 
local spatial correlations are also naturally generated by 
a purely local Hubbard interaction, and, as it is well 
known, they become of essential importance in the vicin- 
ity of second-order phase transitions. As these correla- 
tions are neglected in DMFT, this scheme provides only 
for a conventional mean-field (MF) description of the crit- 
ical properties. 

To overcome this shortcoming cluster extensions to 
DMFT such as the dynamical cluster approximation 
(DCA) and cluster-DMFT have been proposed 0. In 
these approaches spatial correlations beyond DMFT are 
taken into account, however only within the range of 
the cluster size; and due to computational limitations 
the actual size of d = 3-clusters is restricted to about 
100 sites. Hence, short-range correlations are included 
by these approaches, whereas long-range ones are not 
(e.g. for spacings larger than 5 lattice sites in d = 3). 
Nonetheless, Kent et al. 0] were able to extrapolate the 
cluster size of so-called Betts clusters to infinity, albeit 
assuming from the beginning the critical exponents to 
be those of the Heisenberg model. This way they ex- 
trapolated the Neel temperature of the paramagnetic- 
to-antifcrromagnctic phase transition which was found 
in agreement with earlier lattice quantum Monte Carlo 
(QMC) results by Staudt et al. Q. 

As an alternative to cluster extensions and, in partic- 
ular, to include long-range correlations on an equal foot- 
ing, more recently diagrammatic expansions of DMFT 



have been proposed: (i) the DMFT plus spin-fermion 
model [9|, (ii) the dynamical vertex approximation 
(DrA) |10Hl2j which approximates the fully irreducible 
n-particle vertex to be local or that of a DCA clus- 
ter [lH; and (iii) the dual fermion approach [l3j]. As 
for phase transitions, DrA with Moriyasque corrections 
14 1 fulfills - in contrast with dual fermion calculations of 
15| - the Mermin and Wagner theorem in two-dimensions 
and, as we will discuss in the following, corrects the MF 
behavior for the critical exponents in three dimensions. 

In this paper, we apply the aforementioned approxi- 
mation of the DrA scheme (with Moriyasque corrections 
14j ) for studying the phase-diagram of the three dimen- 
sional Hubbard model at half-filling. In particular, we (i) 
calculate the critical exponents, (ii) determine the phase 
diagram with Xjv substantially reduced compared to the 
DMFT one, and (iii) define the region where non-local 
correlations become too strong so that DMFT is not ap- 
plicable anymore. 

We consider the Hubbard model on a cubic lattice 



H = -t 



C i<7 C jv 



(i) 



where t denotes the hopping amplitude between nearest- 
neighbors, U the Coulomb interaction, and c| (cjo-) cre- 
ates (annihilates) an electron with spin a on site i; 
n ia = c \ a c ia ■ I R t nc following, we restrict ourselves to the 
paramagnetic phase with n = 1 electron/site at a finite 
temperature T. For the sake of clarity, and in accordance 
with previous publications, we will define hereafter our 
energies in terms of a typical energy scale D = 2y / 6i flil ]. 

The DTA approach to the model (p]) was derived in 
Refs. [13 and 14. The dynamic non-uniform susceptibility 
reads 
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blc), Gk,i/ = [iv — £k + M — ~£>ioc(v)} 1 (Green func- 
tion) and Sioc(^) (local self-energy). The vertex r^" c " ir 
is determined from the solution of the single-impurity 
problem [To| . In fact, the complete inclusion of non-local 
corrections in the irreducible vertices in all channels can 
be achieved only via the fully self-consistent DrA equa- 
tions. However, as discussed in Ref. 14 1, when con- 
sidering a situation where no competition between dif- 
ferent instabilities occurs, a restriction to one specific 
channel and the evaluation of the self-consistency ef- 
fect via the corresponding Moriyasque correction \ s m 
is possible 14|. In the half-filled case we neglect non- 
local particle-particle fluctuations since this channel is 
strongly suppressed by the repulsive interaction. Fur- 
thermore, at half filling, charge excitations are generi- 
cally expected to be irrelevant for the critical behavior 
as well. Indeed we find Xq W >X{* < x qw (x q £ is the 
particle-particle susceptibility), hence we neglect non- 
local particle-particle contributions as well as A c and de- 
termine A s from the exact sum rule (which also holds for 
DMFT) - ^ImE k>I/ = U 2 n(l - n/2)/2, where the 
non-local self-energy is given by 
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with 7 



s(c),q 



duciblc local spin (charge) vertex, determined from the 
single-impurity problem. 

Starting point of our investigation of the critical prop- 
erties of the antiferromagnetic (AF) instability is the cor- 
responding (divergent) spin susceptibility 



Xaf = Xq,o 



dT(S z , Q (T)S z ,- Q (0)) 



(4) 



with Q = (tt, 7r, 7r). While the DrA with Moriyasque cor- 
rections well reproduces the textbook Mermin and Wag- 
ner results for the Hubbard model in d = 2 yielding finite, 
but exponentially large susceptibility at finite T 14| , the 
situation in d = 3 is even more intriguing, since the AF- 
phase remains stable in a broad region at finite T, allow- 
ing for a direct study of the critical properties. 

Of particular interest is the analysis of the evolution 
of the critical region as a function of the Coulomb repul- 
sion. In Fig. [TJ we show the inverse susceptibility Xaf 
as a function of T for different U values. The vanishing 
of Xaf (T — TV) 7 marks the onset of the AF long- 
range order, defining the corresponding Xjv for a given [7. 
More important is, however, the examination of the crit- 
ical behavior: While in a MF (or DMFT) approach Xaf 
is vanishing linearly close to Tjy in accordance with the 
MF (Gaussian) critical exponent 7 = 1 (see lower inset 
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FIG. 1. (Color online) Inverse AF spin susceptibility as a 
function of T for different U values. Lower inset: Inverse 
DMFT susceptibility with a MF (7 = 1: linear behavior) 
critical exponent. Upper inset: larger T interval. 



of Fig. [1} , DrA data clearly show a bending in the region 
close to the AF transition (i.e., for T < Tq, the so-called 
Ginzburg temperature), indicating a DrA critical expo- 
nent 7 definitely larger than 1. The non-perturbative 
nature of DrA also allows for a treatment of the critical 
behavior, e.g. the size of the critical region, as a func- 
tion of U : From our data it emerges that, in the [/-range 
studied, the size of the region where the critical behavior 
deviates from the MF predictions (here: from linearity) 
increases with U. In order to quantify this statement, 
we have performed DrA calculations at higher T (upper 
inset of FigQ} for U up to 1.5, and fitted the data lin- 
early in the high-T regime. Tq has been hence estimated 
as the temperature below which the relative deviation of 
Xaf fr° m the above-mentioned linear fit becomes larger 
than 10% (red arrows in the upper inset of Fig. [T]). By 
this criterion for Tq-, the size of the critical region with 
non-MF behavior, i.e. AT cr ;t = Tq — TV, increases from 
~ 0.01 for C/ = 1.0, to ~ 0.02 for [7 = 1.25 and ~ 0.025 for 
U = 1.5, following therefore the dependence determined 
by the Ginzburg criterion, which implies the inapplica- 
bility of the standard Landau-Ginzburg expansion in the 
temperature region AT crit oc For [7 < 1 (not 

shown) the bending of Xaf becomes hardly visible, since 
in this regime Tn ~ e~wrj (with W oc 1/D), and there- 
fore the size of the critical region is rather narrow; the 
linear behavior for U> 1.5 becomes confined to temper- 
atures even higher than those shown in Fig. [T] 

A more quantitative study of the critical behavior re- 
quires also a precise evaluation of the critical exponent (s). 
From the behavior of the spin-susceptibility, one can ex- 
tract the values of the critical exponent v, which controls 
the divergence of the AF-correlation length £ (defined as 
the square root of the inverse mass of the spin-spin prop- 
agator at q = Q, u = 0) when T —> Tpj. This can be 
computed cither from the divergence of xaf (i-C-, di- 
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FIG. 2. (Color online) Fit of Xaf( t ) and £ _1 ( T ) for the 
highest interaction value considered, i.e., U = 2.5. Left: fit 
with fixed v = 0.707 (Heisenberg-exponent in d = 3 [2p||). 
Right: free fit, showing the good compatibility with the d — 3 
Heisenberg universality class. 



FIG. 3. (Color online) Neel Temperature in DrA, compared 
with the corresponding DMFT, DCA, QMC and DDMC ones. 
Also shown is the region where corrections to DMFT in form 
of non-local correlations become dominant (violet dashed- 
dotted line: 10 % relative change between the DMFT and 
the DrA self-energies at the lowest Matsubara frequency). 
Inset: DMFT vs. DFA spectral functions (see text). 



rectly from the data shown in Fig. [1} , using the relation 
7 = 2z^ [18( , or by extracting from xaf the value of £ by 
fitting its q-dependence for different T 

The results of our analysis, shown in Fig. [5J demon- 
strate that DrA can describe well the AF criticality of 
the Hubbard model. For the largest values of U = 2.5, in- 
deed, both divergences of xaf and £ observed in DrA can 
be described (left panels of Fig. [2]) with high-accuracy by 
the critical exponent v = 0.707 of the d = 3- Heisenberg 
AF. This is expected to be the correct exponent, not only 
because the half-filled Hubbard can be mapped onto the 
Heisenberg model but also since dimension and symme- 
try of the order parameter suggest the same universality 
class. Similar results, though with a lower degree of pre- 
cision, can be found by directly fitting the value of the 
v exponent to Xaf anc ^ £ (right panels): For U = 2.5, 
our two fits provide an estimate of v ~ 0.70 and 0.73, 
respectively. This shows the Heisenberg universality is 
still valid also in a parameter region (i.e., at interme- 
diate coupling), where the Hubbard model is not well 



approximated by the Heisenberg model [21 



A natural by-product of the calculations of the crit- 
ical exponents is the determination of TV at the DrA 
level, whose values overall well agree with the most ac- 
curate DCA and QMC/DDMC data (sec Fig. EJ). The 
deviations around U = 1 might originate from neglect- 
ing the rather small non-local corrections of the chargc- 
and particle-particle-channels, which could affect non- 
universal quantities such as TV. On the other hand, also 
the DCA/QMC finite-size extrapolation is difficult in this 
region since the AF correlation length is large. Let us 
also note that in this regime the DrA self-energy com- 
pares well with the DCA one of Ref. [HI (U = 1.633, T= 
0.0714): The deviation ± £n |Im E Drv i(k, uj n ) - 



ImE^cACk, w n )|/|ImE,DrA(k,u;„)| is < 5% in the sum 
over the first N = 7 Matsubara frequencies (i.e, for those, 
where a deviation from DMFT is observable). This is 
within the DCA difference between the two largest clus- 
ters considered (84 and 100 sites). 

Finally, we investigate the effects of the non-local cor- 
rections on the spectral properties of the d = 3 Hub- 
bard model. On general grounds, the maximum impact 
of non-local corrections is to be expected close to the 
second-order transition line. This is because the corre- 
sponding spin susceptibility, which explicitly enters in 
the DTA equations for E, is diverging at the transition 
(red line in Fig. [3]) . Such behavior is particularly evident 
in the spectra shown in the two lower insets of Fig. [3] 
for temperatures slightly above the Tn of DrA. Specifi- 
cally, we compared paramagnetic DMFT and DrA spec- 
tral functions at two different k-points on the Fermi Sur- 
face (FS) 0, i.e. ki = (§,§,§), k 2 = (tt,0, §). At 
weak-coupling (U = 1) we observe a strong broadening 
of the DMFT quasiparticle (QP) peak. At [7 = 2, the 
enhanced scattering by non-local spin fluctuations even 
qualitatively changes the spectra: the (already) damped 
QP peak of DMFT is transformed into a "pseudogap" in 
DrA. In principle, one can expect pseudogap behavior 
very close to the Neel temperature also for an arbitrar- 
ily small Coulomb interaction. The corresponding region 
appears, however, at small U very narrow: a qualitative 
estimate according to equ. ([3]) yields the condition for the 
pseudogap behavior £ > 4itvf/(TnU 2 ) (vp is an average 
Fermi velocity), which can be hardly fulfilled at small U, 
where TV is exponentially small. Outside the pseudogap 
region, AF fluctuations yield only an increase of the scat- 
tering rate 7(k) = -ImS(k,w = 0): e.g., at T~Tjy and 
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U = 1 we obtain jdmft = 0.02, jDTA(ki) = 0.033, and 
jDr A {k 2 ) = 0.041. 

By increasing T, non-local corrections become natu- 
rally weaker, since AF-fluctuations are reduced in in- 
tensity and spatial extension, see, e.g., the temper- 
ature behavior of £ in Fig. [2] As a criterion to 
evaluate the impact of non-local correlations, valid for 
the "pseudogap" as well as for insulating spectra, we 
have chosen the relative change between the DrA and 
DMFT self-energy at the lowest Matsubara frequency: 
|£dmft(^i) - S D GA(k 2 , m^)|/|£dmft(m'i)|- Note that 
this criterion is directly related to the QP weight Z in 
the metallic phase if the linear low frequency behavior 
of the self-energy already holds (approximately) at the 
lowest Matsubara frequency iv\. By this one-particle 
criterion, DMFT is reliable down to the violet line in 
Fig. [3] below which deviations exceed 10%. Above this 
line, the impact of the non-local correlations on the 
spectral functions appears indeed moderate (upper in- 
set of Fig. [3]): this is also confirmed by the analysis 
of the spectral function, where the QP weight Z is un- 
changed (within errors) from the DMFT value (Z = 0.76) 
and the enhancement of 7 is much smaller than before 
(idmft = 0.027, lDrA {ki) = 0.028, 7151^2) = 0.036). 

While our findings may validate (a posteriori) the us- 
age of DMFT for computing spectral functions in d = 3, 
provided one is not interested in the immediate vicinity 
of (second-order) magnetic instabilities, it is important 
to note that the width of the critical region is not small 
at intermediate U. For instance, we observe that the size 
of the critical region AT cr u at U > 1.25 exceeds the vio- 
let line. Significant effects of non-local correlations may 
occur even further away from the AF-transition, depend- 
ing on the quantity under consideration. In particular 
relevant deviations from the DMFT predictions at even 
higher-Ts have been reported when analyzing the tem- 
perature dependence of the entropy Q ■ 

In conclusion, we have analyzed non-perturbatively 
the effect of non-local correlations in the d = 3 half- 
filled Hubbard model by means of DrA. When consid- 
ering regions where spatial correlations strongly mod- 
ify the DMFT physics, which is particularly true close 
to magnetic instabilities, DrA represents a very power- 
ful tool for studying the critical properties beyond the 
MF/DMFT level: critical exponents of the Hubbard 
model are found to be -within the error bars- identi- 
cal to those of the d = 3 Heiscnbcrg model, and DrA 
provides also for a proper reduction of TV w.r.t. the 
DMFT prediction. Moreover, since the DrA scheme in- 
cludes both spatial and temporal electronic correlations 
in a non-perturbative way, it looks naturally very promis- 
ing also for future analysis of quantum phase transitions 
beyond the weak-coupling regime. 
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